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Abstract 

We study the exit problem of solutions of the stochastic differential equation dXf = —U' (Xf ) dt+e dLt 
from bounded or unbounded intervals which contain the unique asymptotically stable critical point of 
the deterministic dynamical system Yt — —U'{Yt). The process L is composed of a standard Brownian 
motion and a symmetric a-stable Levy process. Using probabilistic estimates we show that in the small 
noise limit e ^ 0, the exit time of from an interval is an exponentially distributed random variable 
and determine its expected value. Due to the heavy-tail nature of the a-stable component of L, the 
results differ strongly from the well known case in which the deterministic dynamical system undergoes 
purely Gaussian perturbations. 
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Introduction 

The study of dynamical systems subject to small random perturbations keeps receiving much attention 
both in the physical and the mathematical literature. In the simplest (one-dimensional) setting, systems 
of this type find the following mathematical formulation. Consider the ordinary differential equation 
Yt = —U'{Yt), Yo G [—b,a\, a,b > 0, where [/ is a potential function. Assume that U{0) — and that 
is the unique asymptotically stable point of the deterministic dynamical system associated with the 
equation, that means that for any starting point x in [—b,a], the deterministic trajectory tends to 0: 
Yt{x) — > as t — > oo. 

Now perturb the deterministic dynamical system with some small random noise, that is, consider the 
solutions of the stochastic differential equation 

^x- ( U'{Xl)ds + evt, (*) 
Jo 

where 77 is a random process, and the noise intensity parameter e is small compared to the other pa- 
rameters of the system (e is close to 0). Under certain conditions on U and rj, for example under the 
assumptions that U' is Lipschitz and 77 is a semimartingale, the solution of the equation (*) is well de- 
fined. In case ry is a standard Brownian motion, the dynamical system is said to be perturbed by white 
noise. If rj is an Ornstein-Uhlenbeck process, the terminology of 'red noise' perturbation has been used. 
The literature also knows perturbations by the so-called shot noises, fractional Gaussian noises, or Levy 
noises. 

The stochastic dynamics of systems perturbed by white noise, which belong to the large class of 
diffusions driven by the Brownian motion, in the small noise limit, i.e. for e ^ 0, has received a great 
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deal of attention for decades and is particularly well understood. The pioneering papers on this topic are 
IPAV33I |E^r35| IKra40l . Later on it was studied in IFri74l ISctiSOl IWil82l IHsySSl IBM92I IKM96I IFW981 
IKolOOl lBEnTO4l IBGROSl . as weU as in many other papers. 

One of the main results in this field is concerned with the time it takes for the diffusion to exit a 
neighbourhood of a local attractor. Due to the fact that Kramers' pioneering paper was one of the first 
to derive heuristically some properties of an exit law of this type, in particular in the physical literature 
it is often named Kramers' law. Stated in modern language, it says that the expected exit time is 
exponentially large in e^^ and the growth rate can be interpreted as the height of the potential barrier 
to be overcome to leave the local attractor neighbourhood (see section for a rigorous formulation). 

White noise perturbations, however, are not always appropriate to interpret real data in a reasonable 
way. This is the case for example if the nature of the underlying random perturbation process has to 
model abrupt pulses or extreme events. A more natural mathematical framework for these phenomena 
takes into account other than purely Brownian perturbations. In particular infinitely divisible Levy 
perturbations with jumps enter the stage. 

The physical papers by P. Ditlevsen IDit99allDit99bl motivating our research stipulate more general 
noise sources of the type alluded to. They originate in simple physical concepts serving to interpret 
paleoclimatic data. In fact, paleoclimatic records from the Greenland ice-core show that the climate of 
the last glacial period experienced rapid transitions between cold basic glacial periods and several warmer 
interstadials (the so-called Dansgaard-Oeschger events). Those records are given by the concentration of 
certain oxygen, hydrogen or calcium isotopes in the annual layers of the ice-core extending over several 
hundred millennia in the past. They can be used to reconstruct the global Earth temperature for the 
time span for which the records are available. The calcium signal has the highest — almost annual — 
temporal resolution and provides the most conclusive information about the statistics of the Dansgaard- 
Oeschger warming events. They start with a very rapid warming of the North Atlantic region of about 
5-10° C within at most a few decades. The warming is followed by a plateau phase with slow cooling 
extending over several centuries, ending with an equally abrupt drop to basic glacial conditions. The 
nature of these events is not clear, and several conceptual explanations have been proposed. One line 
of arguments invokes instabilities in the North Atlantic thermohaline circulation as a causal mechanism 
for these abrupt climate changes, and for the millennial time scale between jumps. A different reasoning 
starts from the idea that the coupled atmosphere-ocean system in the tropics possesses several meta- 
stable states, and claims global teleconnections that trigger changes in the North Atlantic thermohaline 
circulation. In |GR01I . the effect of stochastic resonance was brought into play in order to explain the 
observed random periodicity of the Dansgaard-Oeschger events. 

As for many climate phenomena, due to the non-linearity of the system the physical background is 
highly complex. To understand some basic features, simple low- dimensional models may be used. For 
the phenomena under consideration, in this spirit one may hope to recover important aspects of the 
statistical properties of the observed data by modelling the paleoclimatic temperature process as the 
solution of an equation of the type (*) with some noise term whose nature has to be determined. To 
take account of meta- or multi-stability, it is natural to assume that the climatic potential function has 
(at least) two wells, their stable minima corresponding to one cold (basic glacial) and at least one warm 
state (plus an intermediate one). In this setting characteristics of the transition mechanism between 
climate states may be reformulated in terms of the exit problem from local attractor neighbourhoods for 
solutions of the stochastic differential equation. This approach was taken in IDit99al |Dit99 bl . To account 
for a reasonable choice of random noise perturbing the system, a spectral analysis of real ice-core data 
was performed in IDit99bl . The obtained spectral decomposition exhibits a strong a-stable component 
with a ~ 1.75. The paper IDit99a| is concerned with an analysis of the exit times of (*) with an a-stable 
noise 77 and in the limit of small e, performed on a physical level of rigour, with the help of a fractional 
Fokker-Planck equation. 

Climate dynamics is not the only source of stochastic models in which a-stable noise appears. For 
example, it was shown in ISMB97I IBS02I . that the thermally activated motion of the test particle along a 
polymer in three-dimensional space is subject to i-stable motion due to the polymer's self-intersections. 
In recent years. Levy noise sources have been playing an increasingly important role in models of financial 
markets (see for example |ER99| 'l. 

In this paper we consider the equation (*) driven by the Levy process L which is the sum of a standard 
Brownian motion and an a-stable Levy motion. Our approach of the asymptotic laws in the small noise 
limit of exit times from bounded intervals or intervals which are unbounded from one side is purely 
probabilistic and completely avoids fractional Fokker-Planck equations. We understand it as a first step 
towards a complete understanding of transition patterns of Levy-driven dynamical systems in bi- or 
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multi-stable potentials. The mathematical challenge consists in a large deviation analysis for exit times 
replacing the classical theory of Freidlin-Wentzell IFW98I for diffusions with Brownian noise. We base it 
on a noise intensity dependent decomposition of L into a sum of two independent processes: a compound 
Poisson with large jumps on the one hand, and a sum of the Brownian motion and a Levy motion with 
small jumps on the other hand. Given such a decomposition for small noise intensity e, the main idea 
of our analysis is to prove that asymptotically exits from the considered domains are due to large jumps 
of the first component, while the second component is not able to perturb the deterministic trajectory 
of solutions of (*) without noise essentially. For this reason, the usual picture of a particle that has to 
climb a potential well being pushed by a Brownian motion in order to exit a domain, which captures the 
system's behaviour for Gaussian noise, changes drastically here. Instead of the height of the potential 
well, a large jump to exit just takes note of the distance from the domain's boundary. Pure horizontal 
distances replace geometric quantities related to the potential in the large deviations' estimates for exit 
times in the a-stable Levy case. Also, the mean values of the exit times change essentially in comparison 
to the Gaussian noise setting: instead of Kramers' times we obtain exit times of the order of in the 
small noise limit, i.e. times of polynomial instead of exponential dependence on e. Our approach can be 
extended to more general heavy-tailed Levy processes. 

The structure of the paper is as follows. In section Q we give a heuristic discussion of the asymp- 
totic exit law based on the decomposition of the a-stable Levy noise perturbing our system, and the 
heuristic picture that the small jump component does not essentially affect the asymptotic behaviour. 
In section 121 we underpin this heuristic picture with mathematical rigour in proving that trajectories of 
the deterministic system and the one in which only the small noise component is admitted to perturb 
are asymptotically very close. This crucial observation is used in sectionj^to derive in a rather technical 
way upper and lower estimates for the law of the exit time of a bounded interval. In sections 3] and |S] 
these results are transferred to the setting of intervals which are bounded on one side. This requires the 
possibility for the deterministic trajectory to return from —oo in finite time. 

Acknowledgements. This work was partially supported by the DFG Research Center Mathematics 
for Key Technologies (Matheon) and the DFG research project Stochastic Dynamics of Climate States. 

Preliminaries and notation 

On a filtered probability space {yi,J-, {J-t)t>o,'P) we consider a stochastic differential equation driven by 
a Levy noise of intensity e: 

X! = x- I U'(XU)ds + eLt, e > 0. (1) 
Jo 

In general, a Levy process is known to be a random process with independent and stationary increments, 
which is continuous in probability and possesses rcU paths. It is completely determined by its one- 
dimensional distributions which are infinitely divisible and characterised by the Levy-Hincin formula. In 
this paper we assume that 

Ee^'^^ =exp|-df + / {e^'y-l^^Xy^\y\<l})-^\, (2) 

[ ^ vk\{o} \y\ ) 

that is L is a sum of a standard Brownian motion with variance d > and an independent a-stable Levy 
motion with < a < 2. More information on Levy processes can be obtained from r5er98.. ,Sat99l . Since 
a Levy process is a semimartingale, the standard theory of stochastic integration applies to equation 
Q, see ISG03llPro04| for more details. Throughout this paper we assume that the underlying filtration 
fulfils the usual conditions in the sense of IPro04l . i.e. the filtration {J-t)t>o consists of cr-algebras which 
are complete with respect to P and is right-continuous. 

The Levy measure of L is given by v{dy) = ^y'^i+c ; y 7^ 0. It is heavy-tailed and has infinite mass for 
all a G (0, 2), due to a strong intensity of small jumps. 

We impose some geometric conditions on the potential function U. First, we assume that U has 
a 'parabolic' shape with its non-degenerate global minimum at the origin, i.e. U'{x)x > 0, U{0) = 0, 
U'{x) = iff 2: = 0, and !7"(0) = M > 0. Further, to guarantee the existence of a strong unique solution 
of Q on R we demand that U' is at least locally Lipschitz and increases faster than a linear function at 
±00 (see also ISG03I fPro04l '). Moreover, in order to obtain some fine small-noise approximations of X"^ 
in section 121 we need that U € in some sufficiently large interval containing the origin. 
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We shall study the first exit problem for the process from bounded and unbounded intervals in 
the small noise limit e ^ 0. In fact, we consider two cases. 

(B) Let / = [—b, a], a,b > 0, and define the first exit time from / as 

aie) = mf{t>0 : ^[-b,a]}. (3) 

(U) Let J = (— oo,a], a > 0, and assume that for some ci,C2 > the regularity condition U{x) = 
Cilxp"'''^^ , X —oo, holds (see Remark 14. II on page 1221 . In this case we study the one-sided counterpart 
of a{e) defined by 

T{e)=mf{t >0 : X! >a}. (4) 

We shall investigate the laws of a{e) and T(e), in particular their mean values, as e — > 0. 

As a main tool of our analysis, we decompose the Levy process L into e-dependent small and large 
jump components. In mathematical terms, we represent the process L at any time t as a sum of two 
independent processes Lt = + Vt , with characteristic functions 



Ee'^«i = exp \-d^+ [ (e'^ - 1 - ^Xy^\y\ < < ^}r;fe 

^ Jm\!(n \y\ 



im\{o} 



(5) 



Ee"^''i=exp / (e-- _ > ^} 

[M{o} ve \y\ 

The Levy measures corresponding of the processes (J^ and 77'^ are 

i.|(-) = K- n{0< |yl< -i=}), ^^(.) = ^(. n{|y|>-L}). (6) 

The process has an infinite Levy measure with support [—^, ■^]\{Q}, and makes infinitely many 

jumps on any time interval of positive length. The absolute value of its jumps does not exceed l/\/e- It 
will be explained later in Remark 13. II on Dage l2Ul whv the threshold is chosen. 

The Levy measure of 77' is finite. Denote 

Then, rj° is a compound Poisson process with intensity and jumps distributed according to the 
law l3-^v^^{-). 

Denote Tk, k > Q, the arrival times of the jumps of r)'^ with tq = 0. Let Tk = Tk — Tk~i denote the 
inter-jump periods, and Wk ~ J?^^ — Vt^- the jump heights of 77'^. Then, the three processes {Tk)k>i, 
{Wk)k>i, and are independent. Moreover, 

P(rfe >u)= I Pee"'^'" ds = e""-", u > 0, 

ETk = -^ = 00 as e -> 0, (8) 



P{Wk £A) = ^ j > -L}^_ dy, for any Borel set A C 

Pe J A \y\ 



Due to the strong Markov property, for any stopping time r the process Qj^-r — Cti ^ ^ 0, is also a Levy 
process with the same law as 

For fc > 1 consider the processes 



x't{x)=x- [ U' {x':_) ds + e^t , t€[0,Tk]. 
Jo 



(9) 



In our notation, for a:: £ R, 

Xt = xl{x) + eWiI{t = Ti}, te[0, Ti], 

= x^A4, + em) + eW2l{t = Ta}, t e [0, T2], ^^^^ 

X!+,^_^ = xUx^;!, + eWk-i) + £WkI{t = Tk}, t G [0, Tk]. 
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Finally, we denote by Y{x) the deterministic function solving the non-perturbed version of Q 

Yt{x)=x~ I U'{Ys{x))ds, x€R. (11) 
Jo 

1 Heuristic derivation of the main result and 
comparison with Gaussian case 

In this section we shall provide the skeleton of a heuristic derivation of our main result on the asymptotic 
law of the exit time from a bounded interval. In the subsequent two sections a rigorous underpinning of 
these arguments will be given. 

On the interval [0,Tfc] and for a; £ R let us consider Y{x) and x''{x). These processes satisfy the 
equations 

Xt{x) = x- U'{x';_{x))ds + e^t, 

t (12) 

Yt{x) = x- f U'{Ys{x))ds, t£[0,Tk], 
Jo 

while the law of Tk is an exponential with parameter /3e . 

The process being the part of L with the small jumps, our analysis will be based on comparisons 
of the trajectories of x'' and Y. If they are close, e.g. if for some 7 > 0, Pc(suPq<3<j,^ jxj — ^s(a;)l > 
is small enough as e — > 0, we can apply the following reasoning. 

For any x £ I, the deterministic solution Y{x) converges exponentially fast to the stable attractor 
0. Define the relaxation time R{x,e) the process Y needs to reach an e'''-neighbourhood of from an 
arbitrary point x £ I. Then, as a separation of variables argument in IIH implies, for some /ii > 0, 

7?(x,e)<max{|J £ < Ine], Q < e < so. (13) 

Since Y{x) and x^{x) are close on the interjump interval [0, Tfe] for all fc > and all x € I, can 
leave I only at one of the time instants while jumping by the distance eW^. 

If the process has not left / at jump number fc — 1, it waits for the next possibility to jump at the end 
of a random exponentially distributed time period Tk . Since Tk (on average) is essentially larger than the 
bound on the 'relaxation' time Ine], e ^ 0, this means that Jf^ jumps from a small neighbourhood 
of the attractor (see Fig.0. 
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Fig. 1: A sample solution of the stochastic difFerential equation 
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Therefore, the foUowing heuristic estimate makes clear, that in the small noise limit cr(e) is an 
exponentially distributed random variable with parameter e°'6 /a, with 



1 1 



(14) 



Indeed, for all k > 1, Tk = X^^^i has a Gamma law with the density p^e (fc-i)! • ^^^'^^ u > 

oo 

P4(j(£) >u)^Y, P(rfc > u) ■ P4o-(e) = Tk) 
fc=i 

oo 

= P(rfc > u) ■ P{eWi £ 7, . . .,eWk-i G /,eVFfc ^ /) 



(fc-l)! 



dt-{l- P(eWi ^ I)Y~'- ■ P(eW^i ^ /) 
(/3,f)'=-^(l-P(eWi^/))'=-^ 



(15) 



/■OO 

J u 

For deriving the last formula we use the equations 



■ dt 



P{eWi ^ J) = P(Wi < or VKi > ^) = 



+ 



1 1 r/e\" e° 



(16) 



The mean value of a{e) may be obtained immediately from 1151 . or independently by the following 
reasoning; 



E,a(e) ^ E^fc • P4«^(e) = ^fc) = I] ^ETfe • P(em £/,..., eVTfe-i G /, eVTfc ^ 7) 



P(eW"i ^ /) 



1 1 



(17) 



The aim of this paper is to make these heuristic arguments rigorous. This is done by proving the 
following Theorems. 

Theorem 1.1 There exist positive constants eo, j, S, and C > such that for < e < eo the following 
asymptotics holds 

exp |-?i£°^(l + C£^)| (1 - Ce^) < P4o-(e) > u) < exp |-?i£"^(l - C£*)| (1 + C£^) (18) 

uniformly for all x G [—6 + £^, a — £^] and u>Q, where = + . Consequently, 

-1 



E,cr(£) 



1 1 



+ )) 



(19) 



uniformly for all x G [—6 + £^, a — £^] . 

Theorem 1.2 There exist positive constants eo, 7, 5, and C > sitc/i that for < £ < £0 t/te following 
asymptotics holds 

exp|-u-^(l + C£*)| (1 - C£*) < P:r(r(£) >u)< exp|-ii-^(l -C£*)| (1 + C£*) (20) 



uniformly for all x < a — e"' and ?i > 0. Consequently, 



E.r(£) = a|^(l + 0(£^)) 



(21) 



uniformly for all x < a — 



It is interesting to compare the results stated above with their well-known counterparts for diffusions 
driven by the Brownian motion of small intensity e. Together with (0 consider the diffusion X"^ which 
solves the stochastic differential equation 



X! ^x- / U'{Xt) ds + eWt, (22) 
Jo 

where W is a standard one-dimensional Brownian motion, and U is the same potential as in (^1. For the 
diffusion X° we define the first exit time of the interval / by 

a(e) = inf{t > : X| ^ [-b, a]}. (23) 

Then the following statements hold for a{e) in the limit of small e. 

1. The first exit time a{£) is exponentially large in e~^. Assume for definiteness, that U{a) < U{—h). 
Then for any 5 > 0,x £ I, according to IFW98I : 

p^(g(2c/w-.)/.^ ^ ,(^)^^(2C7W+6)/.^)_^^ as £^0. (24) 

Moreover, £^lnEa:(T(e) 2U{a). 

The mean of the first exit time can be calculated more explicitly (Kramers' law) |Kra40l ISchSOl 

E.a(e) « £^^_g2C/(.)/.= (25) 

2. The normalised first exit time is exponentially distributed IWil82l |Day83| fbGK05| : for u > 

>v] ^ e"" as £ ^ 0, (26) 



E,a(£) 

uniformly in x on compact subsets of (—6, a). 

As we see, a{e) and a{e) have different orders of growth as £ — > 0. The exit times of the Q-stable 
driven processes are much shorter because of the presence of large jumps which occur with polynomially 
small probability. To leave the interval, the diffusion X^ has to overcome a potential barrier of height 
either U{—b) or U{a). So in the case considered here, = a, with an overwhelming probability. The 

diffusion 'climbs' up in the potential landscape. This also explains why the pre-factor in 12511 depends 
on geometric properties of U such as the slope at the exit point and the curvature at the local minimum, 
the place where the diffusion spends most of its time before exit. 

The process X^ on the contrary uses the possibility to exit the interval at one large jump. This is the 
reason why the asymptotic exit time depends mainly on the distance between the stable point and the 
interval's boundaries. The potential's geometry does not play a big role for the low order approximations 
of the exit time ^{e). Although it is important for the proof, it does not appear in the pre-factors of the 
mean exit time in (1191 and remains hidden in the error terms. 

In the purely Gaussian case, to obtain the law of o-{e), the theory of partial differential equations 
is used. In fact, the probability ps{x,u) = Px{a-{e) > «) as a function of x and u satisfies a backward 
Kolmogorov equation (parabolic partial differential equation) with appropriate boundary conditions. 
The function Ps{x, u) can be (at least in R) expanded in a Fourier series with respect to the eigensystem 
of the diffusion's infinitesimal generator. Then, one concludes that Pe(a;,u) ~ e""^!", where Af is the 
first eigenvalue. Further one shows, that AfE3;(j(£) — > 1 as £ ^ 0. In contrast to this, in the present 
paper we obtain results without any use of operator theory. This suggests that some asymptotic spectral 
properties of the integro-differential operators corresponding to the process X^ can be formulated from 
the probabilistic estimates obtained here. This is the subject of future research. 

Perturbations of the deterministic dynamical systems by small infinitely divisible noises were consid- 
ered e.g. in IFW98| . however in a different setting. There, the small parameter £ was responsible for the 
simultaneous scaling of jump size and jump intensity. As the simplest example of such a perturbation 
one can consider a compensated Poisson process with jump size of the order e and jump intensity of 
the order l/£ (see also IBorGTH . In such a case the dynamics in the limit corresponds to the one of the 
system perturbed by white noise, i.e. the probabilities of the rare (exit) events are exponentially small 
in £ and the characteristic time scales are exponentially large. Note that the perturbations considered in 
the present paper are of quite a different nature. We only scale the jump sizes, while the jump intensities 
stay unchanged. 
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2 Deviations from the deterministic trajectory: exit from 
bounded interval 

In this section we estimate the deviation of the solutions of the stochastic differential equation driven by 
the small-jump process e(,° from the deterministic trajectory on random time intervals of exponentially 
distributed length. We show that the probabilities for at least polynomially small deviations are poly- 
nomially small in e in the small noise limit. This rigorously underpins the starting point of our heuristic 
derivation of the exit law. 

For X £ [—b, a] consider solutions x^ and Y of the equations 

xt{x) = x- U' (xl_{x)) ds + £^t , 

it (27) 



Yt{x) = x- [ U'{Ys{x))ds. 
Jo 



The goal of this section is to prove the following estimate. 

Proposition 2.1 Let be exponentially distributed with parameter (3^, and independent of . Let 
c > 0, 7 = 2-^. Then there is eo > and C > such that for all < £ < £o and x £ [— a] the 
inequality 

P4 sup \xt-Yt{x)\>ce-')<Ce^°'+^'>^' (28) 

holds. 

In order to prove the Proposition, we shall make use of the following Lemma in which the estimation 
of the deviation from the deterministic trajectory is executed on a bounded deterministic time interval. 

Lemma 2.1 Let T > 0, c > and 7 — 2^ (t) < 7 < |J. Then there exist positive numbers eo and C 
such that for all < £ < eq and x £ [—6, a] the inequality 

P^sup \xl - Yt{x)\ > ce^) < CTe"+'^/2 (29) 

[0,T] 

holds. 

Proof of Proposition I^TTl We apply Lemma l2. H and the definition of /3e to obtain for all x £ [— b, a] 

and £ < £0 

r 00 

P^sup |xf-yt(x)| >c£^)= / V^{sxip\xl~Yt{x)\>c£'')p,e-''^^dT 

[O.T.I J, [0..1 

< / rP.e-"'^ dr < C£^^+^'>^\ 



The proof of Lemma l2.1l is performed in three Lemmas in the sequel. It extensively uses the following 
geometric properties of the potential U : 

1. The deterministic trajectories Yt{x), x £ [— &, a] converge to zero as t ^ 00 due to the property 
U'{x)x > for a; / 0. 

2. The curvature of the potential at x = is positive. In small neighbourhoods of zero we have 
U{x) — ^x^ + o(x^). Consequently Y decays there like e^^^*, and the dynamics of x"^ reminds of 
the dynamics of a process of Ornstein-Uhlenbeck type. 

We now prepare our rigorous analysis by an asymptotic expansion of the random trajectories of 2;^ 
around the deterministic one of Y. To this end, fix some 5 > small enough which will be specified later 
and define 

T has the property that for all x £ [—a,b] and t > T, \Yt{x)\ < 5, i.e. after T the trajectory of Y{x) is 
within a (5-neighbourhood of the origin. We next consider the representation of the process x^ in powers 
of £ of the form 

x'{x)=Y{x)+£Z'{x) + R'{x), (32) 
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where is the first approximation of x'^ in powers of e satisfying the stochastic differential equation 

Zl{x) = - / U" {Ys{x))ZU{x) ds + (33) 
Jo 

The solution to this equation is explicitly given by 

Zl{x)= /'e-i"^"(^"(-))*"dC:. (34) 
Jo 

Integration by parts results in the following representation for : 

Zl{x)=et - /*CJ-t^"(n(x))e-J'>"(^"(^»''"ds. (35) 
Jo 

For x = 0, Yt{x) — for all t > 0, and Z^{0) is a process of the Ornstein-Uhlenbeck type starting at zero 
and given by the equation 

ZliO) = Q ~ M f C-e-''^'-'^ ds. (36) 
Jo 

Lemma 2.2 There is a universal constant Cz > 2 such that for any T > 0, x £ [—b, a] and e > 

sup \Zl(x)\ < Cz sup 15," 1, P,-a.s. (37) 

se[o,T] se[o,T] 

Proof: Obviously, for t < T 

\Zt{x)\< sup 15*1(1+ sup /V"(>'.(^))!e" ) . (38) 
te[o,T] y te[o,T]Jo / 

We shall show that the integral in the parenthesis is uniformly bounded. Fix some S > small enough 
such that for some < mi < m2 the inequality mi < ini^^^^^s U"{x) < sup|^|^^- U"{x) < ra2 holds. This 
implies, that mi < U" (Ytix)) < m2 for all x £ I, and t >T. Let 

C7i=max /'^|C/"(y,(a-))|e-^=''^"<^"("»''"ds. (39) 

Jo 

Consider an arbitrary t > T. Then 

/"V"(>'4^))|e" ds 
Jo 

= /"^|f/"(n(x))|e- ''^*^"(^"("»''"ds+ f \U" {Ys{x))\e 
Jo Jf 

Let us estimate the first summand in 1401 . We have 

Jo 



(40) 



, jt u"(Yu{x))du 



ds. 



(41) 



< e-'"i'*-^)Ci < Ci. 



The second summand in l|4()|l is estimated analogously: 



!f/"(y.(x))|e-''='^"(^"(^»''"ds<m2 / e-'"i('-='ds< ^. (42) 

Jf "11 

Taking Cz ~ max{2, Ci + — } completes the proof. 

■ 

To estimate the remainder term R"^ we need finer smoothness properties of the potential U. However, 
the following Lemma shows that this restriction only has to hold locally. 

Lemma 2.3 There exists C > such that for all x £ [—6, a] and T > 0, 

sup \Rl{x)\<C (43) 
te[o,T] 

a.s. on the event {uj : sup^gjQ j,j |e5t('^)i < !}■ 
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Proof: By hypothesis we know that for any t > 0,x £ [~b,a] we have —6 < Yt{x) < a. Moreover, on 
{supjgjQ \sS,t \ < 1} we have supjg[g |e^f | < Cz due to Lemma 12.21 Recall that U' increases at least 
linearly at infinity. This guarantees the existence of C > such that for any y G [—6, a], 2 G [—Cz, Cz] 
we have 

~U'{y + z + C) + U'(y) + U"{y)z < 0. 
Hence for any r>0,T>r>0,a;G [~b, a] the inequality 

-U'(Yr{x) + eZ%_{x) + C) + U'{Yr{x)) + (7"(n(x))eZ;_(x) < 0. (44) 

holds on the event {supjgjQ |£^f| < 1}. Now assume there is some x G [— &, a], and some (smallest) 
T G [0, T\ such that Rt{x) = C . Observe that the rest term satisfies the integral equation 

Rl{x)^ r f{Rl{x),Ys{x),eZ:_{x))ds (45) 
Jo 

with the smooth integrand 

f{R,Y,£Z) = ~U'{Y + eZ + R)+U'{Y)+U"{Y){eZ). 

This implicitly says that i?^ is an absolutely continuous function of time. By definition of r, we have 

< D+Rl{x)\t=r = -U'{Yr{x) + eZ^^_{x) + C) + U'{Y^{x)) + U" {Yr{x))eZ'^_{x) < 0, 

a contradiction. A similar reasoning applies under the assumption R^{x) — —C. This completes the 
proof. ■ 
Lemma 12.31 has a very convenient consequence. It states precisely that the solution process x"^ , 
with initial state confined to [—b,a], stays bounded by a deterministic constant on sets of the form 
{ij : supjgjQ |e^f (aj)| < 1}. Therefore, in the small noise limit, only local properties of U are relevant 
to our analysis. 

We next obtain a finer estimate of the remainder term i?^ on the time interval [O,!']. 
Lemma 2.4 There exits Cf > such that for <T <f 

sup \Rs{x)\ < Cf{ sup leCsD^i Px-a.s., (46) 

se[0,T] ae[0,T] 

on the event {uj : supjQ y] |e^f | < 1} uniformly for x G [—b, a] and e > 0. 

Proof: Using Lemma 12.31 choose K > such that on the event {a; : suptg^ \e^t \ < 1} the processes 
x'^{x),eZ^ (x), R"^ (x) take their values in [—K,K] as long as time runs in [0, T]. For t < T < T, the rest 
term 7?'^ satisfies the following integral equation: 

Rt{x) = /' [-U'{Ys{x) + eZ:_{x) + Rl{x)) + U'{Ys{x)) + U" iYs{x))eZ:_{x)] ds 
Jo 

= - [ [U'{Ys{x)+eZ',^{x)+Rl{x)) - U'{Y,{x)+eZ',_{x))] ds 
Jo 

- f [U'{Ys{x)+eZ;_{x)) - U'{Ys{x)) - U" {Ys{x))eZ;_{x)] ds 
Jo 



(47) 







U"{9l^)Rl{x)ds- / W<-''\e^,_){eZl_{x)f ds 



with appropriate 6l, 6^ £ [—K,K]. Thus 

\Rt{x)\< f L\Rl{x)\ds + \TLCl{sMp\eQ\f. (48) 

Jo [0,T] 

An application of Gronwall's lemma yields the final estimates 

\Rl{x)\ < \TLCle^''{ sup [e^Af < ifLC|/^( sup le^tlf = C^i sup [e^tlf, 
te[o,T] te(o,T] te[o,T] 



sup \Rt{x)\ < Cf{ sup |e^tl)^. 
te[o,T] te[o,T] 

Now we derive a uniform estimate of the rest term R^ on time intervals longer than T. 



(49) 
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Lemma 2.5 There exist positive constants Cr and Ce < 1 such that for T > 

sup \mix)\ < Cr{ sup \eC\f (50) 

se[0,T] sg[0,T] 

on the event 

Et^{lo : sup |e6'l < Ce} (51) 
te[o,T] 

uniformly for x G [— &, a] . 

Proof: Fix some positive T >T and let uj £ {uj : suptg[o |e^f | < 1}. Again, using Lemma l2.3l choose 
K > such that on the event {uj : supjgjgyj |e^f| < 1} the processes x^{x),eZ'^{x), R^{x) take their 
values in [—K,K] as long as time runs in [0,T]. The rest term R"^ satisfies the integral equation 

mix) = R^x) + j' f{Rlix), Ysix),eZl_{x)) ds (52) 

with 

f{R,Y,eZ) = -^U'{Y + sZ + R)+U'{Y)+U"{Y)(eZ). 

Moreover, J?^ is an absolutely continuous function of time. We write the Taylor expansion for the 
integrand / with some 9 £ [~ K,K]: 

f{R, Y, eZ) = ~U'{Y + eZ + R) + U' {Y) + U"{Y){eZ) 

= ~U'{Y) - U"{Y){R + eZ) - ^^^{R + eZf + U' [Y] + U" {Y){eZ) (53) 

^-U"{Y)R-^^^^P-{R + eZf 

Since U G C^, |f/'^>| is bounded, say by L, on [-K, K]. Using the inequality {R + eZf < 2{R^ + e^Z'^) 
we obtain that for all R, Y and Z such that L\eZ\ < A 



f{R,Y,eZ) < -^U" {Y)R + LR^ + L{eZf < -~U"{Y)R + LR' + A = g+{R,Y) 
f{R,Y,eZ) > -U" {Y)R - LR' - L{eZf > -U" {Y)R - LR' - A ^ g~ {R,Y). 

With the help of Lemma 12.21 this immediately implies for T <t < T that 

Rt{x)=R^x)+ [ f{Rl{x),Y,{x),eZl_{x))ds < Rf^ix) + f g+{Rt{x),Ys{x)) ds, 
Jt Jf 

Rl{x)=Rfi,(x)+ f f{Rl{x),Ys{x),eZl_{x))ds>R'f{x)+ f g' {Rl{x),Ys{x)) ds. 



(54) 



(55) 



with A = Z)e^(supjgjg j,] ICt 1)^ a^nd a constant D > 2LC% which will be specified later. 

To estimate the rest term i?^ on [T, T] we apply the subsequent comparison Lemma [2.6l Consider an 
event 

i?i = : \mf{x)\ <^}^{uj ■ Cf{ sup \eet\f < — ( sup \eeAf} (56) 

Thus, setting D > max{2LC|, Cj.m2} we obtain that Pj:(£'i) = 1, 2; G [— b, a] and the conditions of 
Lemma I2. 51 are fulfilled on the event 

E2 = {uj -.ml- 4AL > 0} = : sup < ^^}- (57) 

ie[o,T] 2VLD ^ ' 

On E2 the following inequality holds a.s. for x G [— &, a]: 

on 

sup iJ?j(s)| < — ( sup \eil\f. (58) 
te[f,T] "^2 ts[o,T] 

Thus denoting C/? — max{Cj,, ^^=} and Ce — min{-2^, 1} we may finish the proof. 



11 



Lemma 2.6 (Comparison lemma) Let T >0, Y be a smooth function on [0,T] and Z a roll function 
on [0,T]. Consider the integral equation 

Rt = Ro+ [ f{Rs,Ys,Zs-)ds, te[0,T], (59) 
Jo 

with a smooth function f satisfying 

g~{R,Y) < f{R,Y,Z) <g+(R,Y), R,Y,ZgR, 

g^{R,Y) ^~U"{Y)R±LR^ ±A, R,Y £R, with L,A>Q. 

Moreover, let < rm < U" {Yt) < m2, t G [0,T], and m? - AAL > 0. Then for < t < T the following 
holds: 

1. ifRo<^ thenRt < 
2- 'fRo>~^ thenRt>-2A, 
This yields, that if\Ro\ < then supjgfg \Rt\ < 

Proof: To prove the first statement, together with 1591 consider the Riccati equation 

r+ = Ro + f g+{rt,Y,)ds, t G [0,T]. (61) 
Jo 

Under the conditions of the Lemma, it is enough to prove two statements: 

a) Rt < r+ for < t < T. 

b) r+ < ^ for t > 0. 

To show a) we note that at the starting point i = 0, 

D+Rt\ =[im^:!^-^^f(Ro,Yo,Zo)<g+{Ro,Yo)^rt\ , (62) 
lt=o hio n lt=o 

consequently it follows from the continuity of R and r'^ that r^ > Rt for at least positive and small t. 
Assume there exists r — mi{t > : Rt = r^} such that r G (0,T]. At the point r the left derivative of 
R is necessarily not less than the derivative of r"*" which leads to the following contradiction: 

D-Rt\ - l™ ~ ^ fiR., Yr, Zr-) > r+l 

lt=T hio h \t=T (63) 

/(J?,,n,Z,_) = /(r+,n,Z,„) <,g+(r,,n). 
To prove b) , we compare r+ with the stationary solution of the autonomous Riccati equation 

Pt=Po+ [ {-m2Ps + Lpl + A) ds, t > 0. (64) 
Jo 

Equation 16411 has two positive stationary solutions at which the integrand vanishes: 



Po=P^ = 5(l±jl-^|. (65) 

Repeating the comparison argument used for a), we obtain that if Rq < p~, then rf < , t £ [0,r]. 
Applying the elementary inequalities 



-<l-^/T^<x, a;G[0, 1], (66) 

to p~ we obtain that < p~ < ^ < ■ This guarantees that for Ro < , the solution of l|59|l does 
not exceed — . 

mi 

The proof of the second statement is analogous. ■ 
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Proof of Lemma I2.H Let T > and x G [— &, a]. Choose Et according to Lemma 12.51 Then there 
exists eo such that for e < eo the following holds: 

{ sup \xt{x) -Yt{x)\> ce''} = { sup \eZl {x) + Rt{x)\ > ce''} 
te[o,T] te[o,T] 



C { sup \eZt{x)\ > ^} U { sup \m{x)\ > ^} 
te[o,T] te[o,T] 



C{ sup Htl > #7}U 
te[o,T] ^ 



{ sup \Rt{x)\ >sf}nET 
te[o,T] 



{ sup \Rt{x)\ >sf}nE^r 
te[o,T] 



C { sup Ht] > U { sup |eC?| > fs^} U { sup > Ce} 



te[o,T] 

= { sup Htl > #^}. 

te[o,T] 



te[o,T] 



te[o,T] 



Consequently, with the help of Kolmogorov's inequality we obtain for small e and some C > 



P4 sup \xl - Y,{x)\ > ce-') < P( sup > ^/') < :£|;E(eCT)' 



tS[0,T] 



telo.T] 

,1+q/2-27 .^2-27 



„2^27 



c2 '2 -a 



I+q/2-27 



(67) 



(68) 



This completes the proof. 



3 The law of a{£) 

For the purposes of this rather technical section we introduce the following notation. Denote Wo = To — 

0. a;'^(0) = X, and write lA for the indicator function of a measurable set A. Let denote the exterior 
and interior neighbourhoods of / defined by — [—b =F 5, a ± g] for < q < min{a, b}. 

Throughout this section we use the following constants. Choose fj,i, fi2 and eo > such that for all 
y £ I the following holds for e £ (0, eo], (see Fig.j^J: 

1. |yt(3/)| < for t > /ill lne|. Consequently, denoting z^{y) = x'l:{y) — Yt(y), we have |2;J^(y)| < e"' for 
t > \ lne| under the condition |zj°(j/)| < 

2. Yt{y) £ /j^T for all t > fi2£-' ■ Consequently, |3;t (y)| £ -^3^7/2 foi' * ^ M2e^ under the condition 
[•^t (j/)| < Let us show that /12 with this property exists. Indeed, for all j/ £ /, |yt(y)| is strictly 
decreasing {y / 0), and Yi(0) = 0. Moreover, for — b < yi < 3/2 < a, Yt{yi) < Yt{y2). Let r{x,£) = 
mi{t > : Yt( G ^2£y} • Then by comparison for all x ^ I 

K.,e)<max{y^_^^^^,y_^ ^} < ,2e^ < e < eo. (69) 



3.1 Upper estimate 

In this subsection we give estimates of 'Px{a{e) > u) from above as e — » 0, u > 0. They are comprised in 
the following Proposition with a rather technical proof. Recall that 7 has been chosen according to a in 
Proposition 12.11 

Proposition 3.1 Let 5 = min{a/2, 7/2}. There exist constants eo > and C > such that for all 
< e < eo, X £ [—6 -\- e~' ,a — e''] and u > the following inequality holds 



Pa:(o"(e) > u) < exp < —u — 



1 1 



(l-Ce^) Wl + Ce'). 



(70) 



Proof: For x € /, we use the following obvious inequality 

P^(cr(e) >u)=^ P(rfc > u)Pa;(cr(e) = Tfe) +P^((j(e) > u\a{e) £ (Tk-i,Tk))V^{a(e) £ (rfc_i,rfc)) 
00 

< ^P(rfc > ^*)[p.(a(e) = t^) +V,Me) £ (rfc_i,rfc))]. 



(71) 
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a - 3£T/2 
a - 2e~' 



eT/2 



+ 3eT/2 













































lne| 





















Fig. 2: The dynamics of x^{y) under condition sup \x\{y) — Yt{y)\ < 



Then for any x £ I and G N, applying the independence and law properties of the processes , 
j G N, the following chain of inequalities is deduced which results in a factorisation formula for the 
probability under estimation (compare with CHJ)- 

Px(<7(e) = Tfc) = E,I{A': el,se [0,rfe),X?^ ^ /} 
fc-i 

= n Hxiix'jT'^^ + eWj-i) el,se [0, T,-], 4^ (4-_\ + eWj-i) + eW, G 7} 
< E n sup I{a;i(y) G /, s G [0, T,], ^ (i/) + eW, G /} 



Re 



X sup I{x': {y)(:I,s&[0,Tk],xth, {v) + eWu i 1} 



supl{xi(y) G /,s G [0,Tj],x'^iy) + eWj G /} 



(72) 



X E 



supl{s^(y) G /,s G [0,n],x!i^^(y)+eWk i /} 
.ye-f 



E 



sup \{x\{i)) el,se [0, Ti],x^T, {y) + eWi G /} 



X E 



sup (y) G /, s G [0, Ti] , x't, iv) + eWi i /} 
.ye/ 
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Analogously we estimate the probability to exit between the fc-th arrival times of the compound Poisson 
process 77^, k £ N. Here we distinguish two cases. 

In the first case k = l,x £ Then 

P4a(£) G (ro,ri)) = P4a{s) G (O.Ti)) = E,I{XJ ^ / for some s G (O.Ti)} 



< E 



sup l{xl{y) ^ J for some s G [0,ri]} 



In the second case k > 2,x £ I . Then 

P^(a(e) G (rfc_i,rfc)) = E^I{XJ G /, s G [0,rfc_i],Xj ^ / for some s G (rfe_i,rfe)} 

fc-i 



j=i 



X + £VKfe-i) i I for some s G [0, Tu)} 



< eJJ supl{xi(y) G /,s G {Q,T,\,x'^^{y) + eW^ G /}> 



X suvl{xl-\y) £l,s£ (0,rfc_i],4-_\(y) + £Vyfc_i G /, (y) + eVFfc_i) ^ / for some s G [0,rfc]} 

supl{xj(y) G /,s G [Q,T^],x}r,{y) + eW^ G 7} 
X E supl{xj(y) G /,s G (0,Ti],a;^^(y)+£W^i G /, 2;^(x-:^^ (y) + el^i) ^ / for some s G [0, ^2]} . 

(73) 

Next we specify separately in four steps the further estimation for the four different events appearing in 
the formulae for Pi(cr(£) = Tk) and Pi(cr(e) G (rfc_i,Tfc)). 

Step 1. Consider l{xl{y) £l,s£ [0,Ti\,x\{y) + eWi G /}. For J/ G /, we may estimate 

G /, s G [0, T^],x^T, {v) + eW^ G /} 

= l{xl{y)£l,a£[Q,Tr],x\{y)+eW^£l}\l{ sup > ^} + 1{ sup < ^} 

\ se[o,Ti] se[o,Tii 

<I{ sup \zl{y)\> + sup |^](y)! < G/} 

se[0,Ti] ae[0,Ti] 

= I{ sup \zl{y)\>^} 
se[o,Ti] 

+ I{ sup < ^,:c^^^(y)+eM^i G/}(l{kM/i! < > 

se[o,Tii \ 



<I{ sup > + < Jj} 

se[o,Ti] 

+ I{ sup < ^,|eVKi| > %^,a;^,(y)+£lVi G/,Ti >^iillne|} 

sG[0,Til 

+ ]!{ sup 1^1(^)1 < ^,|£Wi| > ^,a;^,(i/)+£lVi G/,Ti <Ati|lnej} 

sG[0,Til 

< i{ sup |2j(y)| > + HkMl < 

se[o,Ti] 

+ ]I{|£Wi| > sl^eWi G 7+-,}+]I{|eWi| > ^,ri < iii\\ne\} 

= l{ sup > ^f} + I{eVKi G > <7tillnel}. 

sS[0,Ti] 



(74) 
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step 2. Consider l{xl{y) G /, s £ [0,Ti], xip_^{y) + eWi ^ /}. For y G /, we may estimate 
l{xl{y) &I,se [0,n],x'r,{y)+eWi i J} 

<I{ sup > sup \z\{y)\<S^,x\^{y)^eW^m 

se[o,Ti] 3e[o,Ti] 

<I{ sup \z\{y)\ > ^f} 

sS[0,Ti] 

+ I{ sup |4(y)|<^,|eVKi|<^,Ti>/.2e\4,(j/) + em^/} (=0) 

sG[0,Ti] 

+ I{ sup |2](j/)|<^f,|eW^i|<^,Ti<^2£^x:^,(j/) + eVKi^J} 

sG[0,Ti] 
sG[0,Ti] 

+ I{ sup |4(y)|<^,kWi|>%^,Ti<Miilne|,x^,(y)+em^/} 

se[o,Ti] 

< I{ sup \z\{y)\ > ^}+I{Ti < fi2e-'} + I{sWi i J^-, } + j > ^,Tl < fii\\ne\}. 

se[o,Ti] 

Step 3. Consider I{xl{y) ^ / for some s £ [0,ri]}. For y G we may estimate 



^ / for some s G [0,Ti]} < I{ sup \zl{y)\ > ^} + 

sG[0,Ti] 

I{ sup I^Ky)! <^^,a;J(y)^/ for some sG[0,ri]}(=0) 

sG[0,Ti] 



(76) 



Step 4. Consider I{x-^(y) G /, s G {0,Ti], x^^{y) + eWi G /, a;^(a;^^ (y) + eWi) ^ J for some s G [0,r2]}. 
For y G /, we may estimate 

G /,s G (0,ri],4^(y)+eW^i G /, (y) + eW^i) ^ / for some s G [0,T2]} 
= I{a;'(2/) G /, s G (0, Ti],x^t, (y) + eWi G 17-, , (y) + eWi) i I for some s G [0, T2]} 

+ l{xl{y) el,se (0, Ti], x;^^ (y) + eWi G A4"-r , ^U^t, (y) +eWi)iI for some s G [0, r2]} 

< l{x\{y) el,se {0,Ti],x\^^{y)+eWi G 17-,} ■ sup I{xl{y) ^ I for some s G [0,r2]} (77) 

+ G J,s G (0,ri],x^^(y)+£W'i G 

< sup i I for some s G [0, r2]} + (y) G /, s G (0, Ti], (y) + G A4"t}- 

ye 77, 

The first term in the resulting expression in the Step 4 is identical to the expression handled in Step 3, 
while the second term requires an inessential modification of the estimation in Step 1. 

Now we apply H74|l . H75^ . H76|l and H77|l to estimate the expectations treated in Steps 1, 2, 3 and 4 
above. In what follows, c. and C. denote appropriate positive constants. Fix also < (5 < min{7/2, q(1 — 
7),a/2} =min{7/2,a/2}. 

Step 1. Estimate E [sup^g^ I{x-^(y) G /,s G (0,Ti],a;^j + e'PVi G /}] . We get 



E 



supl{xj(y) G /,s G (0,Ti],a;^^ + eWi G /} 



< P( sup \z\{y)\ > + P(£W^i G /+ ) +P(l£lVi| > ^)P(Ti < /ill Inel) 

(78) 



^a/2 



< l_£-/^_(l_Ci£*). 



1 1 

+ 



+ C2£"'^/'-^'lln£ 
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step 2. Estimate E [supy^jl{xl{y) € I,s € {0,Ti],x^^ + eWi ^ I}] . In fact, 



E 



supl{xl{y) e I,s & {0,Ti],x't, + eWi <^ 1} 



<P(sup \zl{y)\ > ^)+P(Ti<fi2e^) 

[O.Til 

+ F{eWi i I--,) + P(leWi| > %^)P(Ti < Inej) 



(79) 



</2 



1 1 

+ 



+ C2e 



,(3/2-7) 



llnel 



Step 3. Estimate E j^sup^^^- \{x\{y) ^ J for some s e [0,Ti]}j . We have 



E 



sup I{xl{y) ^ / for some s £ [0,Ti]} 



< P( sup \zl{y)\ >^)< Cie("+^)/2 < e"''^ ■ C,e' (80) 

[0,Ti] 



Step 4. Estimate 

E [sup^g^ ll{a;s(y) G /, s G (0, Ti], 2:^^ (i/) + eWi G /, a;^(a;i,^ (y) + eWi) ^ I for some s G [0, T2]}] . We fi- 
nally obtain 



E 



supll{xj(y) el,se {0,Ti],XT^{y)+eWi G I , x^,{xT^{y) + eWi) ^ I for some s G [0,T2]} 
ye/ 



<P(sup 14 (y) I > %^) 

[O.Ti] 

+ P(sup \zl{y)\ > ^) + P(£Wi G [-fe-£^,-fe + 2e'']) + P(eWi G [a-2e^a + £^]) 

[0,Ti] 

+ P(|eWi| > %^)P(Ti < Mil Inel) C4e'. 
Then for x £ JJ^ , < e < eo, and some positive C5, 



(81) 



P:,(cr(e) > u) < P(ri > u) 



+ 5]P(r, >u)(l-e'^''2^(l-Ci£^)) e' 



1/2 1 



l-ea/2|(i_C'^£i) 



+ 1 + 026" 



{k-l)\ 



2 

1 + Cse^ 
1 - Cse^ 



exp < — ue™ — (1 — Cse ) > < exp < — ite™ — (1 — Cse ) (1 + Cee 



< exp <; -ue"^(l - Ce") [. (1 + Ce") 

(82) 

with C = max{C5,C6}. 

In the previous formula we have changed summation and integration. This can be done due to 
the uniform convergence of the series J2T=i fte~'^'* ''fSiy^ [l - £"''^|(1 - C'se'')] for all t > and 
e < eo- Indeed, let tl be the coordinate of the maximum of the density of the (/?£, fc)-Gamma distribution. 

k-l 



For fc > 2, it is easy to see that tl 

< f3,e 
Then, for all e < eo, 



(k-iy. 



Then, with help of Stirling's formula we obtain, that 

fc-i 



^^^^„(.,i) (fc-l) 



-At 



{k-iy. 



l_e-/^l(l_C5e ) 



{k-iy. V2^Vk^' 



< Cl- 



k 00. 



< c, 



V2|(l_C5£^) - ' 



(83) 



(84) 
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where the constant c does not depend on t and e. Thus uniform convergence foUows from dominated 
convergence. 



3.2 Lower estimate 

In tliis subsection we estimate Px{(j{e) > it) from below as e ^ 0, u > 0. This leads to the following 
Proposition with a rather technical proof again. 

Proposition 3.2 There exist constants eo > and C > such that for all < e < eo, < 5 < 
min{a/2, 7/2}, x G [— fe + e^, a — e~'] and u> the estimate 



(l + Ce')} {1-Ce') 



e° 1 1 
is valid. 

Proof: We use the following inequality: 

oo 

P^(cr(e) > u) = ^P(rfe > u)P^(cr(e) = r^) +P^(cr(e) > u\a € (rfe_i, rfe))P^((j £ (rfc_i,rfc)) 
fc=i 

oo 

>^P(rfc >it)P4a(e) = rfc). 
1 

With arguments analogous to l|72|l we obtain the factorisation 

= Tk) = E^I{XI el,se [0, r,-), x^, ^ /} 
fe-1 

= J] II{2;^,(a;^ir_\ + eW^j-O G s G [0, Tj], a:^^^. (4-_\ + eWj-i) + eW^ G /} 



> I{xi(4-_\ + eW^,_i) G /, s G [0, T,], 4 (4"^ + eW^,_i) + eWj G /^^ 

X I{4(4:_\ +eW/fc-i) G /,s G [0,Tfc],4jx^;_\ +£W'fc_i) +£Wfc ^ /} 

fc-i 

> E n inf_ Hxi {y)eI,se{0,T,], 4 (y ) + eW, G J,". } 



Re 



X inf_ l{x':{y) el,se (0, T,.], (y) + eWk i /} 



inf l{xl{y) el,se (0,T,],4.(j/) +eW^,- G I^-,} 



X E 



inf I{4(y) G /,s G (0,r,],4,(y) ^ /} 



E 



inf G /, s G (0, Ti], 4i (y) + eWi G I^^ 



X E 



inf (y ) G /, s G (0, Ti] , a;^^ (y) + eW^ i 1} 



(85) 



(86) 



(87) 



For y G I^,-,, we next specify separately in two steps the further estimation for the two different events 
appearing in the formulae for Pa;(cr(e) — Tk). 
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(88) 



Step 1. Consider the event I{xl{y) £ I,s £ (0, Ti], (y) + sWi G I^-y}- We may estimate 
l{xl{y) £l,s£ (0,n],x'r,{y)+eWi £ 1^.} 

>I{ sup \zl{y)\<sl,xl(y)eI,se{0,Ti],x^T^{y)+£WieI--,} 
se[o,Tii 

> I{ sup \zl{y)\ < ^,xl{y) el,se iO,Ti],\eWi\ < ^,Ti > /iae^ x^^^ (y) + eW^i G 4".} 
se[o,Tii 

+ I{ sup \zl{y)\ <sl,xl[y)el,se {Q,T{\,\eW^\ > ^.Ti > I Ine], x;^, (y) + eW^i £ 

sG[0,Ti] 

>I{ sup \zl{y)\<sl,\£Wi\<sl,T^>^2e'} 

sG[(),Ti] 

+ I{ sup |4(y)|<^,|£VKi|>^,Ti>Millnei,eVKieJ2-,} 

ae[O.Ti] 

= II{i£Wi| < %^,Ti > ^l2e''} 
-l{ sup I^Ky)! > sl,\eWi\ <^,Ti> fi2£^} 

ae[O.Ti] 

+ I{leWi| > ^,Ti > iii\lne\,eWi G I2,-,} 
~l{ sup 1^1(^)1 > JJ,|eW/i|>^,Ti> Mil lne|, em Gla".} 

sG[0,Ti] 

^I{\£Wi\ < sl}-l{\£Wi\ < ^,Ti<fi2£^}-I{ sup \zl{y)\ > ^,\£Wi\ < %^,Ti > fi2e''} 

se[o,Ti] 

+ I{|em| > -^,eWi G ^"t} -HleW^il > %^,Ti < Mi|lne|,£Wi G Ij^} 

-I{ sup |4(y)| > %^,|em| > Y'^i > Millnei,eM G /2>} 
se[o,Ti] 

>I{£Wi G/a'.^l-IlTi <M2£^}-2I{ sup \zl{y)\ > si} -i{\eWi\ > Ti < Inej}. 

sg[0,Ti] 

Step 2. The event I{xl{y) G /, s G (0, Ti], a:;^^ (y) + eWi ^ 1} may be estimated as follows: 
I{xl{y) el,se (0, Ti], x^, (y) + eW^i ^ J} 

>I{ sup \zl{y)\<sl,xl{y)&I,se{0,Ti],x'TAy)+<^Wi<^I} 

ae[0,Til 

>I{ sup \zl{y)\<^,xl(y)eI,se{0,T-i],Ti>^ii\\ne\,x\,^{y) + eWi<^I} 
se[o,Tii 

>I{ sup \zl(y)\<^,n>^^i\lne\,eWi(^I+} 
3e[o,Ti] 

= I{ri >Ati|ln£|,eWi ^ /+}-!{ sup > ^^.Ti > ^il lne|,eWi ^ /+ } 

se[o,Tii 

> I{e Wi ^ J+ } - I{ri < ^ti 1 In e| , £ W'l ^ /+ } - I{ sup ] ) ] > ^ , eWi i it-, } . 

se[o,Ti] 

Now we apply (1881 and 1891 1 to estimate the expectations appearing in the formula for Pa;(cr = t^). In 
what follows, Ci and d denote appropriate positive constants. Fix also < 5 < min{7/2, 0(1 — 7), = 
min{7/2,a/2}. 

Step 1. Here we estimate E j^inf^^^- I{xl (y) G /, s G (0, Ti], + eWi G I^j }j . In fact, we obtain from 
employing results from sections |5| and 



(89) 



E 



inf I{xl (y) G /, s G (0, Ti] , x\-^ + eWi G } 



> FieWi G 12.-,) - P(Ti < M2e^) - 2P( sup \zl{y)\ > ^) - P(|eWi| > < Mil lne|) 

[0,Ti] 



(90) 



> 1 - 



1 



+ 



1 



(a-2£^)° (fe-2£^)" 



Cl£"/^+^ - 2C2£("+^'''^ - C3£"''/'-^'|ln£| 



> 1_£"/2Z(1 + C'l£*). 
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step 2. We next estimate E [inf^g,- I{xl{y) el,se {0,Tj],x^^ + eWi i /}] , for which we obtain 
similarly 



E 



inf \{x\{y) £l,s£ (0,Ti],4^ +£Wi ^ /} 



PieWi ^ J+ ) I 1 - P(ri < Aill lne|) - P( sup \zl{y)\ > 

lO,Ti] 



(91) 



i/2 t 



{a + e-')°' (b + £T)^ 



1 - C2£ 



1 — cse"''^! In el — C2£ 



(q+7)/2 



Consequently, with C3 = max{Ci, C2}, for a; £ , < e < eo 



1/2 
1 - Cae 



(fc-1) 



l_e-/2^(l + C3e') 



e ' -(1 - Cse ) 



> 



^ exp<^-«e'^-(l + C3e')^ >exp<^-ue"-(l + C3e')[-(l-C5e') 



(92) 



> exp<^ -ue"-(l + Ce*) } {l~Ce 



with C — max{C3, C5}. See the end of the proof of Proposition 13.11 for the justification of switching the 
order of summation and integration in the above argument. 

■ 

Proof of Theorem ll.il The statement of Theorem ll.il follows directly from Propositions 13 . 1 1 and 13.21 



P{eWi i \-h,a\) 



(93) 



Remark 3.1 The threshold l/\/e for separating the two parts of L is not the only possible choice. 
Indeed, let us consider thresholds of the type 1/e'' for p > 0. Then = > j^}u{dy) = ^e"''. p 

must satisfy some further conditions. Firstly, we demand that p < 1 so that we can easily calculate the 
probability 

This is the characteristic probability of our analysis. The probabilities of other relevant events should 
have smaller order in the small noise limit. This for example applies to the event that leaves the 
e'^-tube around the deterministic trajectory (see Proposition 12. H . and the event that Ti < fj,2£'') (e.g. 
see 1781 . 1791 1 and 1I8UI1 ). This leads to the following inequalities on p and 7: < p < 1, 7 > 0, 
a{l — p) < 2 — 2p — 27 and a(l — p) < ap + 7. Applying some algebra we rewrite these inequalities in 
the form 

\<2^{l-p), 0<p<l, 
7>a(l-2p), 7 > 0. 

The solution set (p, 7) is non-empty for all a G (0, 2) and depends on a. However p = i is the minimal 
value independent of a for which there exist 7 solving the inequalities. In this case any 7 from the 
interval (0, is a solution. For our purposes we have taken 7 — ^-k^- 



(94) 



4 Return from — oo and deviations from the deterministic 
trajectory: exit from unbounded interval 

With the aim of proving an analogue of Proposition 12. II we study in this section the exit problem of the 
solution of our stochastic differential equation from the unbounded interval J = (— 00, a]. In this case we 
shall use the condition that U' increases faster than a linear function at —00 which guarantees a return 
from infinity in finite time for the unperturbed deterministic motion. For simplicity we assume that for 
large |a::|, [/ is a power function, i.e. |?7(a;)| = ci |a;|^''"'^^ , ci,C2 > for \x\ > N. This condition can be 
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weakened (see Remark 14. II at the end of this section). We stick to it to avoid technicaUties irrelevant for 
the main aim of the paper. 

Fix two more positive numbers r and R such that N < r < R, and such that with Tr = iw^y)] 
we have — r — YTj^i—R). Consider the equations 1271 on the unbounded interval J. As in section |5| 
we estimate on the basis of the representations x'^{x) — Yt{x) + eZt{x) + Rf{x), with Yt{x) = x — 
U'{Ys{x)) ds and Zt{x) = ~ /g' £,3-U" {Ys{x)) ^/|y*(!^j) ds. We first prove an estimate enabling us to 
transfer Lemma 12.21 to unbounded intervals. 

Lemma 4.1 The inequality sup^gjQ j,^] |Zt (a::)| < 2sup(gjQ | holds a.s. for x < —R. 

Proof: For all x < —R by definition of R and r we have Yt(x) < — r. Moreover, by assumption 
U"{Yt{x)) > for t e [0,Ta], whence 

Zax)< sup 16^1(1+ sup sup r U"{Y,{x))%^}^^ ds] . (95) 
te[o,THi y ^<-Rteio,Tn]Jo ^ ('^3{^)) ) 

We show that the integral in the latter parenthesis is uniformly bounded in x. Denote Yt{x) — v. Then 
dv = Yt{x) dt = -U'{Yt{x)) dt = -U'{v) dt. Therefore 



= U'{Y{x)) 



,2 
1 



U'{Yt{x)) U'{x) 



< 1. 



Here is the analogue of Lemma 12.21 
Lemma 4.2 There is a constant C'z > such that the inequality 

sup \Z't{x)\<C'z sup m (96) 
te[o,T] te[o,T] 

holds a.s. for x < a, T > 0,e > 0. 

Proof: The proof obviously has to combine the previous Lemma with Lemma l2.2l on bounded intervals. 
The inequality holds with C'z = 1 + Cz ■ B 
To estimate deviations of the random paths from the paths of the deterministic equation, we restrict 
from the start to sets of bounded scaled noise. 

Lemma 4.3 On the event {sup^g^ je^f | < 1} the following inequality holds a.s. 

sup \xt - Yt{x)\ < 10 (97) 

uniformly for x < —R. 

Proof: It follows from Lemma l4.ll that supjg[Q \eZl{x)\ < 2 a.s. on the event {supjgjg |e^f | < 1}. 
We show that the rest term \R°\ is bounded by 8. Indeed, the rest term satisfies the integral equation 

RUx) = I [~U'{Ys{x)+eZl^{x) + RUx)) + U'iYs{x)) + U"iY,{x))eZl_ix)] ds. (98) 
Jo 

R^{x) is absolutely continuous a.s. and Ro{x) — 0. Assume, that there exists a smallest r G [0,Tfl] such 
that R%{x) = 8. Then the left Dini derivative of R^{x) at this point is non-negative, i.e. 

-U'{Y^{x) + eZ'^_{x) + 8) + U'{Yr{x)) + U"{Yr{x))eZ',_{x) > 0. (99) 

On the other hand, our conditions on U guarantee 

- U'{Y^{x) + eZ%^{x) + 8) + U'iyrix)) + U" {Yr(x))eZ%_{x) 



< ^U'{Yr{x) + 6) + U'{Yr{x)) + 2U"{Yr{x)) < 0, 



(100) 



and a contradiction is reached. The estimate Rt(x) > —8 is obtained analogously. ■ 
The following Lemma is an analogue of Lemma f2 . 3 1 and gives a rough estimate for the remainder term 
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Lemma 4.4 There exists C > such that for x G [—R, a], T > 

sup \Rl\ < C (101) 
te[o.T] 

a.s. on the event {supjgjg j,j ]e^f| < 1}. 

Lemma f4.4l again has localising consequences: It states precisely that the solution process x"^ , with ini- 
tial state confined to [—R, a], stays bounded by a deterministic constant on sets of the form {sup^gjQ |e^f (lj)| < 
1}. Therefore, in the small noise limit, only local properties of U are relevant to our analysis. 

Let us paraphrase the most important aspect of what we found in the previous Lemmas due to finite 
return from — oo. With probability close to one, the random trajectory starting ai x < —R reaches the 
finite interval [—R, —r] in a finite non-random time Tr which does not depend on e. Our investigation 
therefore reduces to the study of the dynamics of paths starting in the finite interval [—R,a]. Since 
the deterministic trajectories starting in [—R,a] do not leave this interval, the statement of Lemma 12.21 
does not change. Due to Lemma 14.41 the estimate of the rest term R"^ given in Lemma 12.51 also holds 
unchanged. 

Thus, in the case of the unbounded interval we have all necessary tools to estimate the exit proba- 
bilities. 

Remark 4.1 The conditions on the behaviour of the potential U at — oo can be weakened. Indeed, a 
slight extension of the proof of Lemma k. II allows to drop the convexity condition U" > 0. Furthermore, 
to show that 7?^ is bounded from above by some pi > 0, we need to guarantee that the integrand in 
(I98II is negative for R^ — pi under the condition \eZ^^\ < p2, P2 > 0. This leads to the inequality 
— infeg[p2_pj_p2+pi] U'{y + 6) + U'{y) +pi\U"(y)\ < 0, which has to hold for y < —N. For instance, for 
the power function U{x) — cilip^'^^ considered above this inequality is equivalent to — (p2 — 2pi)|y|~^ + 
C(|2/|~^) < 0, y —00, and thus holds for any p2 > 2pi > 0. 



5 The law of t{£) 

In this section we estimate Pi(r(e) > u) for w > as e — > 0. Indeed, the extensions of Propositions 13. ll 
for the estimation above and 13.21 for lower bounds with parameter b — +oo take the following form. 

Proposition 5.1 Let 5 = min{a/2, 7/2}. There exists eo > and C > such that for all < e < eo, 
x £ (—00, a — e^] and u > 

P4r{e) >u)< exp|-u-^(l-C£*)| (1 + Ce*). (102) 

Proposition 5.2 Let S < min{Q/2, 7/2}. There exists eo > and C > such that for all < e < eo, 

X G (—00, a — e^] and u > 

P.(r(e) > u) < exp|-u-^(l + C£*)| (1 - Ce^). (103) 

Proof: The arguments to prove the estimates are similar to those of the bounded case. We just need 
to adapt Steps 1, 2, 3 and 4 (see sections 13.11 and 13.2^ to the case of an unbounded interval. Let us 
consider for example the extension of Step 1 from section I^TTI The basic formula 1721 holds in the case of 
unbounded intervals with J replacing I. We demonstrate how to modify the reasoning just in the series 
of inequalities 1741 and in the estimate (I78II . The other estimates are obtained analogously. 

Firstly, we estimate I{xl{y) G J, s G [0,Ti],x\,^{y) + eWi G J} for -R < y < a. Denote A = 
{suptg[Q j,^] \e^t \ < !}• On the event A, the trajectory xl{y), t G [0, Ti], belongs to a compact interval, 
so its dynamics is indistinguishable form the one treated in the bounded case. Therefore, we have 

l{xl{y) G J, s G [0, Ti],x^T^ (v) + eWi £ J} 
= I{xliy) G ^,^6 [0,T,],x'T,{y)+eWi G J} {I{A} + 1{A^}) 

< I{xl (y) G J, s G [0, n],x}r, (y) + eWi eJ,A} + I{A'} <■■■ (^O^) 
<I{ sup \xl{y)-Ys{y)\ > sl} +i{eWi G J+ } -f I{|eW^i| > ^,Ti < ^,i\\ne\} + I{A''} 

[0,Ti] 
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The case y < —R is slightly more complicated since we have to treat the return of x^{y) to a compact 
interval in a finite time. Denote 



B = {uj : sup |e^t I < f} n {i^ : sup |eCt - £^ThaTi1 < 1} 

[O.Tj^ATi] [ThATi.Ti] 

D {iu : sup \eet \ < 1} ^ A. 

[0,Ti] 

Then we have 

I{xl{y) e J, s G [0, Ti],x't, iv) + eWi G J} 

< Hxliy) e J,se [0,Tr ATi],xl{y) £j,se [Tr AT^,T^],x\-^{y) + eWi G J}nB + I{A"} 

< l{xl{y) eJ,se [Tr a TuT^],x^r, (y) + eWi G J} n B + I{A''} <■■ 

< I{ sup \xl{y) - n_T„ATi(a^T^ATi(y))l > + I{eM^i G J+-,} 

se[THATi,Ti] 

+ I{|eVFi| > ^^.TflATi < ^ijlne|}+I{A"} 

< sup I{ sup \xliy) - Ys{y)\ > ^} + l{eWi G J+.} 

+ I{|eM| > ^,Ti < ^i|lne|} + I{A'=}. 
These estimates may be treated in a way similar to H78|l . In fact, 



(105) 



(106) 



E 



supl{a::j(j/) G J, s G (0,Ti],a:^^ + eWi G J} 



< 1 



7^(1 -Ce") 
2a" ^ ' 



(107) 



for some C > 0. The other steps are modified analogously. 

Proof of Theorem 11.21 Combine the estimates of the above Propositions. 
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